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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMDRANDUM 

THEORETICAL ANALYSIS OF THE LONGITUDINAL BEHAVIOR OF AN 
AUTOMATICALLY CONTROLLED SUPERSONIC INTERCEPTOR 
DURING THE ATTACK PHASE AGAINST MANEUVERING 
AND NONMANEUVERING TARGETS 
By C. H. Woodling and Ordway B. Gates, Jr. 

SUMMARY 


A theoretical analysis has been made of the longitudinal behavior 
of an automatically controlled supersonic interceptor during the attack • 
phase of the interception problem. Attack runs were computed for a 
nonmaneuvering target, and for a target which had a constant acceleration 
normal to its flight path. First-order lead collision navigation was 
ass\jmed in the investigation, and characteristics of this navigation when 
used against a maneuvering target are discussed. The flight path of the 
interceptor was controlled by commanding either a pitching velocity or 
normal acceleration proportional to the vertical steering error. Com- 
puted attack runs are presented which demonstrate some of the advantages 
and disadvantages of using high gain or integration in the tracking system 
to minimize or eliminate bias errors in the system which result from tar- 
get acceleration or interceptor trim changes. 

Results are also presented which show the effect of limits on the 
rate of control deflection, and several means of co-unteracting the effects 
of this limiting are discussed. 


INTRODUCTION 


The Langley Laboratory of the National Advisory Committee for 
Aeronautics is conducting, an analytical study of automatically controlled, 
rocket firing, supersonic interceptors during the attack phase of the 
interception problem. The term "attack phase" refers to that phase of 
the interception which exists subsequent to the time at which the inter- 
ceptor's radar becomes locked on to a specified target. Results have 
been reported in reference 1 for the case in which the. interceptor locked 
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onto a nonmaneuvering target with an initial vertical tracking error. 

Only longitudinal maneuvers of the interceptor were necessary to carry 
out the desired interception. Lead collision navigation was used in 
this reference^ and the tie-in between the radar -guidance con 5 >uter and 
the interceptor was a command on pitching velocity proportional to the 
existing error in the interceptor's flight path. The interceptor con- 
sidered in reference 1 has .a notched delta wing of aspect ratio 5.2 
with 55° sweepback of the leading edge. A more detailed discussion of 
this configxiration can be found in reference 2. 

The present paper is concerned primarily with the longitudinal attack 
performance of the interceptor of reference 1 against a target which is 
assmed to be maneuvering with a constant normal acceleration. The inter- 
ceptor is initially in level flight at 50,000 feet at a Mach number of 2.2. 
The target is initially in level flight at a Mach nvimber of 1.4, flying 
toward the interceptor at various altitudes above 50,000 feet. No con- 
sideration is given to the effect of altitude changes on the interceptor's 
longitudinal behavior during the attack runs. 

The results of this investigation are presented, for the most part, 
in the form of interceptor and kinematic responses subsequent to radar 
lock-on, which were computed on the Reeves Electronic Analog Computer 
(REAC) . 


SYMBOLS 

ly moment of inertia about Y-axis, slug-ft^ 

m mass of airplane,, slugs 

c mean aerodynamic chord, ft 

S wing area, sq ft 

q dynamic pressure, Ib/sq ft 

V forward velocity, ft /sec 

M Mach number 

n normal acceleration, number of g 

An change in normal acceleration, number of g 

g acceleration due to gravity, ft/sec^ 
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q 0e/2V when used as a subscript 

9 angle of pitch, radians unless otherwise specified 

a angle of attack, radians unless otherwise specified 

7 flight-path angle (7 =0 - a), radians unless otherwise 

specified 

u or AV change in forward velocity, ft /sec 


u' 
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e 
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T 

e 
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L 


C 


L 


D 


"D 




C 


m 


relative change in forward velocity, H or ^ 
elevator deflection, radians imless otherwise specified 

time, sec 
lift, Ih 

trim lift coefficient, 

^ qS 

drag, Ih 

trim drag coefficient, 

qS 

pitching moment, Ih-ft 

pitching -moment coefficient, Pitching moment 

qSc 







aa ' 


per radian 


per radian 


per radian 


per radian 
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D 

a 


R 

ft 


T 


M 


M. 


LS 


^NLS 




per radian 


— per radian 

a ®£ 

2 V 


differential operator, ^ implies integrationj 

angle between interceptor X body axis add radar line of sight, 
positive when line of sight is above body axis, radians 
vinless otherwise specified 

distance from interceptor to target along line of sight, meas- 
ured positive from interceptor to target, ft 

angvilar velocity of , line of sight, (ft = ct + 0 ), radians/sec; 

. positive when line of sight is rotating upward 

time of flight of interceptor from instantaneous position to 
firing point, sec 

time of flight of interceptor's rockets from firing point to 
predicted point of contact with target, sec 

predicted miss distance, measured positive from interceptor 
to target, ft 

component of M along the instantaneous radar line of sight, 
positive when target is ahead of rockets at predicted time 
of impact, ft 

component of M perpendicular to the instantaneous radar line 
of sight, positive when target is below rockets at predicted 
time of impact, ft 

error in interceptor's flight path at any given instant, 

£ « ~^NLS 

^ VjtQ + + VrJ t 

y 

elevator servo-system time constant, sec 


€f output of filter, radians 
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ratio of steering-error differentiation gain to steering 

. ^ 

error gain — , sec 
K 

Subscripts : 

I interceptor 

R rocket 

T target 

L limit 

ss steady state 

i input 

o initial value 

A dot over a quantity indicates differentiation with respect to time. 
A bar over a quantity indicates a vector. 
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DISCUSSION OF GUIDANCE EQUATIONS, LONGITUDINAL CONTROL 
SYSTEM, AND INTERCEPTOR EQUATIONS OF MOTION 


Guidance equations . - The longitudinal tracking problem against a 
target maneuvering with a constant normal acceleration Vip7,ji is shown 

diagrammatically in figure 1 for a lead collision type of navigation. 

For this type of navigation the interceptor endeavors to fly^a straight- 
line path such that at only one point on the path the rockets may be 
fired and a hit obtained on the target. The vector equation relating 
the distances travelled by interceptor, rockets, and target in the time 
interval (t(j + x), subsequent to any given instant, to the predicted 
miss distance is; 

R + 2 sin 7ip ^ ^ ^ ^ ~ ^I^G ^ ^ ^ (l) 


where R is the range vector from the interceptor, which is flying with 
the velocity Vj, to the target, which is flying with the velocity V,p. 

The parameter is the time rate of change of the target flight path, 

and Vrji/7ip is the radius of curvatxire of the target flight path. The 
time of flight from the given instant to the point at which the rockets 
are released is tg, and x is the time of flight of the rockets from 
the firing point to the predicted point of contact of the rockets with 
the target. For this investigation x = I.5 seconds. The average rocket 
velocity, due to its own thrust, over the time x is Vp. The distance 
between the rockets and target at the predicted time of impact is the 
miss distance vector M. The components of equation (l) along and per- 
pendicular to the instantaneous line of sight are: 


R + 2 


^ sin 7 t cos [a + 0 - ( . ^G + t)J ^ 

IJi^G (^I Vfj^'^cos(a + a) + 

2 sin 7m '^) sin [a 

7T ^ L 

[^1 % + (^I + ^r)"^ sin(a + a) + Mj^g 


|o + e - 7 t - ■j'T 


( 2 ) 
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For the assumption that 


cos 7n 


fa t) ^ 

r\ 


and 


sin 


7t 


(^G + ^ fa + 

. 2 2 


these equations become: 


R + V^^tg + t) cos (a + e - 7,p) + I sin^ci + 0 

^I (% '^) cc.s(a + a) + Vj^T cos(a + a) + M^g 





V,j,(tQ + t) sin^a + 6 - 7^) - J \7^, (tg + cos + 6 - 7^J = 
^l(^G + a) + VpT sin(o + a) + 


(3) 


This assxmiption essentially means that the components of the target 
normal acceleration along and perpendicular to the instantaneous radar 
line of sight are near-ly constant over the time internal (tg + t) . Equa- 
tions (3) rewritten in terms of range rate R and line of sight rota- 
tion n are 


R + R^tG + t) - Vj^T cos(a + a) + i V,p7,j, ^tg + sin + 0 - 7,pj = M^g 
-Rfi(tg + x) - Vpx sin(q + a) - i (t^ + x)^ cos (a + 0 - 7^) = Mj^g 


where 


R = Vip cos + 0 - 7^^ - Vj cos(ct + a) 
-RG = Vpi sin^q + 0 - 7p^ - Vj sin(q + a) 


The target flight -path angle 7^ is taken as zero when the target is in 

level flight going away from the interceptor, and as « when in level 
flight coming toward the interceptor. For 7p < «/2, a positive value 
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of indicates a pitch-up of the target flight path, and for 7,^ > n/2, 

a negative value of 7^ indicates a pitch-up. For the rtins in this paper 

the initial flight-path angle of the target is taken as it (target coming 
toward interceptor) . 

In reference 1, the tie-in between the guidance and the longitudinal 
control system was accomplished by con^juting continuously the value of 
(tfj + t) necessary for M^g to be zero,, and for these values of (t(j + t), 

computing the predicted value of Mjjj^g. The command to the control system 
was based on the error which exists at any time in the interceptor's 

flight path, and is approximated by 

■' 1 % + ("l + Vr)t ^5) 

The time at which tQ is computed to be zero is taken as the firing time 
for the interceptor's rockets. Examination of equations (4) indicates 
that it is necessary to know the target flight-path angle 7tj and target 

normal acceleration Vip7,p in order to solve these eqmtions for (tg + t) 
and MjjLg. In order to evaluate the acceleration terms of equations (4) 
in a guidance computer, knowledge of derivatives of R of higher order 

than R is required. Guidance systems of the type from which these 
derivatives would be available are difficult to mechanize, and are not 
being used. In view of this condition the target acceleration terms in 
equations (4) will be omitted, and the primary effect^ of the maneuvering 
target is assimied to be reflected in the parameters R and Rft through 
changes in the target flight-path angle 7iji. When the acceleration 

terms are omitted, equations (4) become the first-order guidance equa- 
tions presented in reference 1 for a nonmaneuvering target. It is appar- 
ent from examination of equations (4) that the Vip7ip term caji have only 

a negligible effect on the computed value of (tQ + t), and can result 
in a maximvim error of 56 feet in Vqi 7^ = 52 ft/sec^ and 

T =1.5 seconds, which are the values assimied in this investigation. 

One of the most significant characteristics of this first order 
guidance is that the interceptor must maintain, in the steady state, a 
normal acceleration proportional to that of the target in order to carry 
out the tracking assignment against a target maneuvering with a constant 
normal acceleration. For the pitch-rate command system used for longi- 
tudinal control in reference 1, an error in the interceptor flight -path 
angle is required to command a steady interceptor normal acceleration, 
and hence this system is unable to track a maneuvering target with a 
zero error. The miss distance associated with this required 
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flight-path error will henceforth be referred to as the bias error. The 
introduction of a slow integration in the tracking loop should permit 
tracking with a zero bias error. Another possibility is to use a high 
static sensitivity between the interceptor normal acceleration and the 
interceptor flight -path error angle, which should minimize the bias error. 
However, both conditions should be investigated with respect to stability 
and performance against nonmaneuvering and maneuvering targets. 


Control systems considered .- The automatic control systems considered 
in this investigation for controlling the interceptor’s flight path include 
the pitch-rate system discussed in reference 1, and a normal acceleration 
control system which will be discussed subsequently in this paper. Block 
diagrams for both systems are presented in fig\ire 2. The dynamics of the 
filter and elevator servo are represented by first-order lag networks of 


the form 


1 + TfD 


and 


■= respectively. For this investigation 

1 + TgD 


Tf = 0.60 second and Tg =0.05 second. The filter would be used in 
practice to filter the noise out of the error signal No considera- 


tion is given to noise in this investigation, but the assimied dynamics 
of the filter are included in the analysis. The error signal sub- 

sequent to being filtered, is passed through an amplifier, a differen- 
tiator, and an integrator. The result of these operations is taken as 
the command to either a pitch-rate control system or normal -acceleration 
system. In either system the interceptor normal acceleration may be 
limited by limiting the input to these systems to the desired value. 

For this investigation the inputs to the control system were limited to 
values such that the interceptor static acceleration response would not 
exceed +5g or -2g. 


Interceptor equations of motion .- The longitudinal dynamics of the 
interceptor were represented in this investigation by the following equa- 
tions of motion, referred to wind axes: 


— 7 = /Ct Ai + C, be\(l + u') + 2Cj u' + — (A6 - £a) 

W cos 7 ^ 

I . ' O/ 


W sin 7 q 


mv j 


-(A0 - Aa) 


( 6 ) 




q^Sc 


“ (•^q ® ^ (1 + 2u' ) 


% 


From mpublished wind-tunnel tests -made for a model similar to the inter- 
ceptor discussed in this paper, the variation of drag coefficient Cj^ in 


CONFIDENTIAL 



10 


CONFIDENTIAL 


NACA RM L55G18 


the vicinity of the interceptor ' s trim angle of attack and initial Mach 
number was found to be well approximated by 

^D \ 

= 0.027 + 0.156 Zki + 2.57 - (0.015 + 0.151^ At + 2.05 Ax2)u' 

. and the variation of with Mach number in this range was given by 

P _ ^- 0 ^ _ 2.29 

Mq(1 + u') 1 + u' 

The interceptor stability derivatives, mass characteristics, and other 
constants used in this investigation are presented in table I. A detailed 
derivation of equations (6) is presented in appendix A, and the assimqjtions 
made are discussed. 


RESULTS AND DISCUSSION 


In the section entitled "Discussion of Guidance Equations, Longi- 
tudinal Control System, and Interceptor Eq\iations of Motion" it was stated 
that an Interceptor, which utilizes first-order guidance, must have a 
steady normal acceleration proportional to that of the target in order to . 
track a target maneuvering with constant normal acceleration. Furthermore, 
for the longitudinal control systems discussed in this paper there must 
exist an error in the Interceptor ' s flight path, with respect to that 
required for a hit, of sufficient magnitude to command this acceleration 
unless there is integration performed on the flight -path error. The 
magnitude of this bias error can be minimized by use of a high sensitiv- 
ity between the interceptor's g-response and flight-path error if no 
integration is included. It may be necessary, however, to make additional 
modifications to the automatic-control system or tracking loop in order 
to insure adequate stability when high gains are used. 


Control Systems 

Pitch-rate command system . - The block diagram of the pitch-rate 
command system is shown in figure 2(a) . In reference 1 good tracking 
performance against a nonmaneuvering target was calculated for this sys- 
tem for K = 5-0 = 0.575- For a maneuvering tanget, however, 

it can be shown that there will be a large bias error in the flight path, 
and hence miss distance, for this value of K. Since (Zhi)„„ = — 9„oi 

SS g S S ’ 
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the following expression, which relates the interceptor steady normal 
acceleration to the value of Mj^g at the time of firing (tg = 0), may 

he derived from the block diagram of figure 2 and equation (5); 



tG=0 


-gT 

KV-r (1 + u') 
o 




+ Kj. 


^ss 


( 7 ) 


For 

T = 1.5 

Vj + Vp = 4,140 

^ = 0.015 
Kj. = 0.375 

K 3 = 1 


(based on two degrees of freedom) 


(An)ss « An^ = 1 


to=0 ' ^ 

Therefore, for Mjj^g to be less than 50 feet, when tg = 0 , for this 
assvraied case K must be approximately 10 or greater. For K = 5.O, which 
was used in reference. 1 for the nonmaneuvering target, g would be 
approximately -I60 feet for this target acceleration. For a practical 
case in which the interceptor speed would decrease, gains even larger 
would be required. 



and if u' = 0 
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System response without integration in tracking loop: Results are 

presented in figure 5 which afford a comparison of the interceptor's 
attack performance against a nonmaneuvering target and against one maneu- 
vering with Ig normal acceleration when no integrator is in the tracking 
loop. The r\ms shown in figure 5(a) are for control system gains which 
give a reasonable response for the case of a nonmaneuvering target. For 
these runs Rq = 60,000 feet, = 7-5°^ 7j = 0, and 7 ij^ = jt. For 

these and subsequent runs the transient responses are plotted up to the 
assumed time of rocket firing (tg = 0). For comparison, lines of constant 
flight -path error of 20 mils are shown. These lines are shown only in 
this figiire but are the same for all the rvins presented. It is apparent 
that for = 1 , the predicted value of is about -500 feet when 

the rockets are fired (tg = O) . The value of at t(j = 0 predicted 

by equation ( 7 ) for (An)gg and u'gg shown on figxrre 5 (a) is about 

-260 feet, as compared to the value of -500 feet indicated by the REAC 
calculations. There is a small bias error attributable to the intercep- 
tor's loss in forward speed which can be seen from the results for the 
nonmaneuvering case. As the interceptor loses speed, a continuous error 
is generated in the flight -path angle and, probably more important, a 
change in the trim angle of attack occiors as the interceptor speed is 
reduced. This bias can be eliminated, or at least markedly reduced, by 
introduction of a signal to the elevator servo proportional to the loss 
in forward speed, which corrects for the out-of-trim pitching moment due 
to the loss in speed. Results for the case where this feedback is added 
are also presented in figure 5(s-) • This bias due to u' could also be 
eliminated by use of an integration between and 6 g, or minimized 

by use of a high gain between these quantities . The responses presented 
in figure 5 (b) are for a. set of control-system gains which were about the 
best found for the maneuvering target case. The responses for both tar- 
get conditions are seen to be very oscillatory for this high gain and no 
improvement was obtained from further increases in K^. 

System responses with integration in tracking loop: Results are 

presented in figures k(a) and 4 (b) which show the effect of integration 
in the tracking loop on the system responses against maneuvering and 
nonmaneuvering targets . The effect of the integrator is to eliminate the 
bias error against the maneuvering target (fig. 4 (a)), but the integration 
causes a large overshoot in the miss distance for the nonmaneuvering tar- 
get case (fig. 4(b)). Another disadvantage of the integrator can be seen 
from figiire 5* Responses are shown for various initial lock -on errors, 
and although the responses are good for 0 ^ = 7 * 5 °> “tbe integral gain 

used is too large for Oq = 10° and too small for © 0 = 2 °. It appears 
that if an integrator is used it may be necessary to use an integrator 
gain which is a function of lock-pn error, or some similar nonlinear 
arrangement, in order to get satisfactory responses over the range of 
likely to be encountered. It should be pointed out that, when a high 
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gain is used between 0^ and instead of the integration, the quality 

of the system responses are essentially independent of Oq. 

Effect of high gain plus differentiation in tracking loop: The high 

gain cases presented in figure 3(b) were seen to be very oscillatory and 
hence undesirable. This condition covild be improved by the use of dif- 
ferentiation in the tracking loop to provide flight-path stabilization. 

The equation for would then become 


0j^ = Ke^ + 

Resiilts are presented in figvires 6(a) and 6(b) for maneuvering and non- 
maneuvering targets which show the effect of the differentiation gain K2 
on the system responses. There is a large increase in tracking loop 
damping as Kp is introduced. Although K2 tends to stabilize the 
tracking loop, it tends to reduce the effectiveness of the filter, and 
hence might not be acceptable from the standpoint of system noise. This 

©i 

can be seen by considering the transfer function — which is 

^ ^ K(1 + TdP) 

(1 + xfD) 

where 



Therefore, the introduction of differentiation in the tracking loop is 
somewhat comparable to reducing the filter time constant. For exanple, 
when = Tf, the filtering is completely eliminated. It can be deduced 

from these resiilts that this high-gain system could be made very stable 
provided that a filter time constant roughly 50 percent of the value 
ass\mied (0.60) woiild give acceptable filtering. 

Normal -acceleration command system . - In the present investigation, 
the basic error involved is e^, which is the instantaneoxis error in the 
interceptor's flight path. The interceptor normal acceleration is pro- 
portional to the rate of change of the flight path; hence, a system in 
which a normal acceleration is commanded proportional to the flight-path 
error appears to be a very logical type of longitudinal control system 
for the present problem. Control systems similar to this system are 
discussed in references 5 and 1 *:. The block diagram of the normal- 
acceleration system is presented in figure 2 . It should be noted that 
there is an integration in the normal -acceleration command loop. A 
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signal equivalent to cos 0 is subtracted from the accelerometer output 
in order to zero the feedback signal to the elevator servo when no change 
is desired in the flight path of the interceptor. Subtraction of this 
quantity effectively converts the accelerometer output to a An or 
feedback. 

System response without integration in tracking loop (K^ = 0) : 
Resxilts are presented in figure 7 for nonmaneuvering and maneuvering 
target conditions for two values of steering-error gain K. The results 
shown in figure 7(a) are for K = 55 -0^ which gives a good response for 
the nonmaneuvering target, but there is seen to be a large bias error in 
the miss distance for this gain for the maneuvering target case. It will 
be noted that there is no bias due to u' for the nonmaneuvering case’ 
for this control system such as that computed for the pitch-rate system. 
This bias is eliminated by the self -trimming properties of the normal- 
acceleration system which are provided by the integrator in the normal- 
acceleration command loop. The responses shown in figure 7(1^) are for 
K = 200.0. The responses for both target conditions are oscillatory but 
the bias error for the maneuvering target case is reduced to a fairly 
low level at the assuned firing point. For this high value of K it 
was necessary to reduce the gain of the integrator in the control loop 
to zero in order to keep the overall system stable. This removes the 
self -trimming properties of the system, but the bias error due to change 
in trim is very small since the static gain between and &g is 

= 42.0 for this case. 

System responses with integrator in tracking loop: The effects of 

integration in the tracking loop on the system responses are seen in 
figures 8(a) and 8(b) . The same trends are noted for this normal- 
acceleration system as were noted for the pitch-rate system. The 
responses for the maneuvering target case for = 7*5° with K = 35*0 
and = 0.95 (fig* 8(a)) indicate that the bias error is eliminated by 
the integration, but for the nonmaneuvering target condition (fig. 8(b)) 
there is a large overshoot in the transient. Although results 

are not presented for values of Oq other than 7-5°^ 'the same dependence 
of the motions on lock-on error exists for this system as for the pitch- 
rate system discussed previously. 

System responses with high gain plus lead in tracking loop; Results 
are presented in figures 9(a) and 9(b) which show the effect on the sys- 
tem stability of incorporating lead or differentiation in the tracking 
loop, and the results shown for this control system are similar to those 
obtained for the pitch-rate system. The differentiation in the tracking 
loop has a large stabilizing effect on the system responses, but as 
mentioned previously, it might not be desirable from the standpoint of 
system noise. 
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Effect of Limiting Rate of Control Deflection 

The responses presented up to this point were computed for the 

assumption that (6=) = ±120°/sec and = ±20°. In order to 

' '=^max ' ' ®^max 

investigate the effects of reducing the maximum available control rate, 
calculations were made for (^e) may “ 80 °/sec and 60 °/sec and the results 

are presented in figures 10(a) and 10(b) for two values of pitch-rate gain 
Kj.. The results are for the normal-acceleration control system, and are 
for high gain only in the tracking loop. For Kj. = 1, the system is seen 
to be unstable for = 6o°/sec. However, for K^. = 2.0 (fig. 10(b)), 

this instability is not present. This comparison indicates that pitch- 
rate feedback cein be lised to eliminate instability caused by low aAmil- 
able control rates. The results presented in figvire 10(c) demonstrate 
another means of eliminating instability caused by low available control 
fates. For this run, (Be)max “ 6o°/sec and a feedback to the servo 

proportional to pitching acceleration has been assumed. This feedback 
is seen to have a very stabilizing effect on the condition shown. The 
runs presented in figure 10 for Kj. = 2.0 were repeated for the case 
where the steering-error differentiation was included in the tracking 
loop, and the results are shown in figures 10(d) and 10(e) for K 2 = 120. 

It can be seen from figure 10(d) that inclusion of lead resvilted in insta- 
bility for (5e)niax = 6o°/sec, whereas for = 0 (fig. 10(b)), the 
system was stable for this control rate. However, inclusion of pitch 
acceleration feedback (fig. 10(e)) is seen to eliminate this instability. 


CONCLUSIONS 


From the results presented in this paper the following conclusions 
can be drawn: 

1. For the control systems considered, a bias error will exist in 
the flight path of an interceptor which utilizes a first-order guidance 
computer that must accelerate in the steady state when tracking a target 
which is maneuvering with a constant normal acceleration, unless an 
integration is performed on the flight-path error. 

2. The effect of integration in the tracking loop is to cause large 
overshoots against nonmaneuvering targets and to give a system response 
against maneuvering targets which depends on the initial lock-on error 
and target accelerations. 

5* Use of high tracking gain in lieu of integration tended to mini- 
mize the bias errors against maneuvering targets but the stability of 
the system responses, particularly those obtained when the pitch-rate 
command was utilized, was poor. 
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U. Inclusion of a signal proportional to the derivative of the 
filtered flight -path error as part of the command. to the control system 
resTilted in good tracking stability for both the pitch-rate and normal- 
acceleration systems. This type of signal, however, might be undesirable 
from other considerations such as system noise. 

5 . Bias errors in the flight path which result from an interceptor 
trim change dviring the attack rvm, can be eliminated by use of integra- 
tion in the control loop, or minimized by use of a high gain between the 
flight -path error and the elevator deflection. Also, a feedback to the 
elevator servo proportional to change in forward speed was shown to be 
capable of eliminating or minimizing the bias errors due to trim changes. 

6. Reductions in the available control-deflection rate were shown 
to have a destabilizing effect on the system responses, but this desta- 
bilizing effect, could be counteracted by \ise of high pitch-rate feed- 
back or by use of pitch-acceleration feedback. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 

Langley Field, Va., July I 5 , 1955- 
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APPENDIX A 


DERIVATION OF AIRFRAME LONGITUDINAL EQUATIONS USED IN INVESTIGATION 


The airframe equations of longitudinal motion used in this paper 
are derived from the following general longitudinal equations which are 
referred to wind axes: 


mV 7 = L - W cos 7 
mV = -D - W sin 7 + T 
lyG = M' 


(Al) 


The thrust T is assumed to he aligned with the wind at all times. 
Substitution of 


Ij = Lq + 

7 = 7q + A7 


T = T. 


V = V^ + AV = Vq( 1 + u') 

D = Dq + AD 
\ M' = 

L^ - W cos 7=0 

U O 

-D^ - W sin 7 q + Tq = 0 

into equations (ai) yields the following equations, for the assumption 
that cos A 7 = 1 and sin Ay = A 7 : 
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mV^Cl + u')7 = AL + A7W sin 7^ 
mVj^u' = -AD - A7W cos 7^ 


(A2) 


lyG = m' ^ 

The quantities AL, /S), and AM' may be obtained from the relations 


(Lo + AL) = + ^l)(<1o + ^)S 

(Do + ad) = + ACj)^ + Z!>q.)s 


AM' = ACjjj^q^ + Aq)Sc 
The expressions for these quantities are 


AL = Cl Z!qS + ACLqoS + 


AD = C]-, + Z^pq^S + ACj) Z^S 

o 


AM' = ^mqoSc + AC^ AqSc 


(A3) 


The quantities C^ and Cjj are functions of the trim angle of attack 
aQ and the initial forward velocity Vq, and q^' is a function of Vq 

and the initial altitude of the interceptor. The quantity Aq is given 
by 

Ai = qo(2u' + p' + 2p'u' + u'^ + p'u'^) (aL) 


where ' 


and 



u' 


= ^ 
"Vo 
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Substitution of equations (AJ) into equations (A2) yields the following 
set of equations: 


mVf, , 
-^(1 + 

u')r = 

Cl 

^L 

/ 

(1 + 




Co qp 


\ 

%J 


' - c 

— + Z5Cts 



w 

4_ 

‘loS 

1^0 

<lo 

r 

%) 

T — 

ly .. 

= |1 





q^Sc 

m 1 
V 






(A5) 


For the assumption of no change in density, equation {Ak) becomes 

— = 2u' + u'^ 




Substitution of this expression into equations (A5) gives 




-I 


Iq r |2u' + u'^) W sin 


raVr u 


-^ g — = Cp (2u’ + u'2) + ^j3(l + 2u’ + u'2) + A7 

Ho^ o qQO 


— ^ e = zc„(l + 2u' + u'2) 

qQSc 

which, when only first-order terms in u' are retained and for 


a 


'V= 


2v_ ® 2V ^ ^ ^ ^ 

^ In S 
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reduced to 


mVi 


q S 


/ ^ ^ sin 7- 

L. 1 L_ L,, el' V n5(l+u') 


mV 


I0S 


In • W COS 7_ 

° u' = 2CnU' + ^J)+ — — „ ■ — ■ A7 






V 


*q 2Vj 


which are the equations presented as equations (6) of the report. 
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TABLE I. - STABILITY DERIVATIVES AND MASS CHARACTERISTICS OF 
INTERCEPTOR AND OTHER CONSTANTS USED IN INVESTIGATION 


Altitude, ft . . 

p, slvigs/ft^ . . 

Vt , ft /sec . . 
■^o 

m, slxigs .... 

p 

ly, slugs -ft 

q, Ib/ft^ . . . 

c, ft 

S, sq ft . . . . 
Cj„ , per radian 

Cj^, per radian 

C™., per radian 
a 

Cjj ^ , ) radian 

Ci^, per radian 

Cj) , per radian 
a 

"1-0 

Cl 

C~ , per radian 
% 


Cig , per radian 


Vr, 

ft /sec 

Vt, 

ft/sec 

^s^ 

sec 

Tf, 

sec 

T, 

sec . . 


ft . . 

CD 

0 

V* 

radians 


radians 


Kg, radians/radian/sec 



. . . . 50,000 

. . 0.0003622 
2 l 40 (M = 2.2) 

. . . . 776.4 

. . 2.68 X lo 5 


826 

15 

401 

-2. a 


-0.56 


-0.28 


0.00 


2.29 

0.156 

0.027 

0.076 

-0.295 

0.165 

2,000 

1,360 

0.03 

0.60 

1.5 

60,000 

0.033 

0.033 

1.0 

4.85 
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Figure 2.- Block diagrams of the longitudinal tracking systems used in 

present investigation. 
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(b) Normal -acceleration cQiamand system. 













interceptor attack performance with integration in tracking loop 
Pitch-rate command system; Rq = 60,000 ft; K = 3*0; = 0.575j 
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Figure 7 .- Interceptor and kinematic time histories. Normal-acceleration 
command system; Rq = 60,000 ft; Qq = 7 . 5 °; 7 -r = 0 ; 7 m = jt. 
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Figure J.- Concluded. 






Figiire 8.- Concluded- 
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Figure 10.- Continued. 
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Figure 10.- Concluded 





